QCD inequalities for the nucleon mass and the free energy of baryonic matter 
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The positivity of the integrand of certain Euclidean space functional integrals for two flavor QCD 
with degenerate quark masses implies that the free energy per unit volume for QCD with a baryon 
chemical potential /i_g (and zero isospin chemical potential) is greater than the free energy with 
isospin chemical potential ^/ = 2ia. (and zero baryon chemical potential). The same result applies 
to QCD with any number of heavy flavors in addition to the two light flavors so long as the chemical 
potential is understood as applying to the light quark contributions to the baryon number. This 
relation implies a bound on the nucleon mass: there exists a particle X in QCD (presumably the 
pion) such that Mjv > ^Tx^ where mx is the mass and of the particle and Ix is its isospin. 
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PACS numbers: 12.38.Lg,14.20.Dh,21.65.+f 

Quantum chromodynamics (QCD) is the theory under- 
lying strong interactions. The theory is not analytically 
tractable via perturbative and other weak coupling meth- 
ods except for a limited set of observables in a limited 
kinematic regime. While certain aspects of the nonper- 
turbative regime of QCD have been explored via lattice 
QCD.lJ, it is important to establish as many properties 
of QCD as possible via rigorous analytic means. QCD 
inequalities represent an ideal method to do this. An 
early variant of the approach was based on a demon- 
stration by Nussinov 2] that bounds could be placed on 
hadronic quantities in a wide class of models inspired 
by QCD. The approach was greatly strengthened by the 
realization of Weingarten ,3j and Witten 4] that simi- 
lar bounds could be obtained directly from QCD itself 
through inequalities based on Eulcidean space functional 
integral representations of physical quantities. Thus, the 
approach is a method to deduce certain qualitative fea- 
tures of QCD from first principles. While the resulting 
inequalities have not been proved as theorems at the level 
of rigor demanded by mathematicians, they make use of 
only the most standard assumptions made by physicists. 
For example one assumes that the theory exists, that 
physical quantities may be represented via functional in- 
tegrals, that a Wick rotation from Minkowski space to 
Euclidean space is innocuous and so on. The field of 
QCD inequalities is now two decades old; the state of 
the art is presented in a recent comprehensive review by 
Nussinov and Lampert 9. Ahhough the resuhs of QCD 
inequalities are qualitative, they can serve to supplement 
understanding gleaned from lattice QCD studies; they 
provide an analytic means of understanding some of the 
results of QCD which are both seen in nature and which 
emerge from numerical studies on the lattice. QCD in- 
equalities may also provide insight into certain proper- 
ties of QCD which are not tractable on the lattice using 
Monte Carlo algorithms. 

The underlying idea of QCD inequalities is quite sim- 
ple. One relates a physical quantity of interest to a Eu- 
clidean space functional integral over gauge field config- 



urations. If a second quantity can be represented via 
a functional integral for which the integrand is greater 
than or equal to the integrand of the initial quantity for 
all gauge configurations, then one can conclude that the 
second functional integral is bigger than the first and this 
in turn allows one to bound one physical quantity by an- 
other. 

This letter focuses on the use of QCD inequalities in 
two apparently unrelated problems in strong interaction 
physics. The first is bounding the mass of the nucleon 
from below in terms of other physical observables. This 
is an old problem. Nussinov derived a bound that the 
nucleon mass must be greater than or equal to ^^^^ in 
the context of QCD-inspired models. Weingarten, in his 
original paper, attempted to bound the mass of the nu- 
cleon as some multiple of the pion mass directly from 
QCD j^. Unfortunately, the only rigorous method found 
required the study of QCD with six or more degenerate 
light fiavors. An alternative approach was also suggested 
by Weingarten; however, it used plausible but not prov- 
able assumptions about bounds on the quark propagator 
in the presence of an arbitrary gauge potential. Sub- 
sequently, Nussinov and Sathiapalan showed that in 
the large Nc limit of QCD with two degenerate fiavors, 
the nucleon mass is bounded by Mjq > ^"^^ . However, 
previously no rigorous lower bound has been obtained for 
the nucleon mass directly from QCD at finite N^. for two 
degenerate flavors. Such a bound is derived here. 

The second issue discussed in this letter is the problem 
of finding a lower bound on the free energy of QCD at 
nonzero baryon chemical potential, or more precisely a 
baryon chemical potential associated with the two light 
degenerate flavors. This second problem is significant be- 
cause for the zero temperature case it is in the class of 
problems for which Monte Carlo algorithms cannot be 
used in lattice simulations. Moreover no viable alterna- 
tive presently exists for doing such simulations. Thus any 
reliable results from the theory are extremely welcome. 
The problem is also interesting in light of the intense 
recent interest in QCD at finite baryon density 0- 
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The bound on the QCD free energy at fixed baryon 
chemical potential will play an essential role in bound- 
ing the nucleon mass and accordingly this problem will 
be treated first. This problem is rather unusual for a 
QCD inequality treatment. The approach is more com- 
monly associated with bounding masses of particles via 
the study of correlation functions rather than with ther- 
modynamically intensive quantities such as free energy 
densities. However, QCD inequalities have been used 
for intensive quantities in the past. Vafa and Witten 
demonstrated that the vacuum energy for QCD including 
a 9 term has an absolute minimum at = 0. While the 
validity of a related argument by Vafa and Wittenj^ that 
parity cannot be spontaneously broken has recently been 
questioned , the demonstration that the the minimum 
vacuum energy is at = is clearly correct. 

The Vafa- Witten proof is very simple. The integrand 
of the functional integral is given by ni=flavors det(^ 
-I- rui) e-'5i'M +»«^ where det(^-|- rrii), the functional de- 
terminant for a given flavor is real and non-negative 0, 
Sy m is the Yang-Mills action, and v is the winding num- 
ber. Thus, the only effect of setting 9 to be nonzero is 
to multiply the rest of the integrand (which is real and 
positive) by a pure phase factor e^^" . Since the real part 
of this phase factor is always less than or equal to unity 
(and the imaginary part will integrate to zero) one im- 
mediate deduces that the functional integral for nonzero 
9 is bounded from above by the integral with 9 = 0. Fi- 
nally, identifying the functional integral as the generating 
function Z{9) = e^^^^^*^ where V is the four dimensional 
volume and £{9) is the vacuum energy as a function of 9, 
one sees that the inequality for the generating function 
implies that £{9) > £{0). 

Here an argument analogous to that of Vafa and Wit- 
ten will be given for the problem of two flavor QCD with 
degenerate quark masses at a nonzero chemical poten- 
tial. The free energy density for QCD at fixed temper- 
ature and baryon chemical potential, Gb{T, ij-b) is given 
in terms of the grand partition function Zb{T, (Ib) as 
Qb{T,iib) = - if3V3)-^log{ZB{T,^lB)) where is the 
three dimensional volume of the system and P is the in- 
verse temperature. In QCD with two degenerate flavors, 
Zb{T, /is) can be represented as a functional integral. 
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where Nc is the number of colors (3 for the physical 
world), the functional determinant is over a single quark 
flavor and the temperature is implemented by imposing 
periodic boundary conditions in time for the gluon fields 
A{t + (3) = A{t) with (3 = 1/T; similarly antiperiodic 
boundary conditions for the fermions are imposed in the 
functional determinant. The factor of simply reflects 
the fact that the chemical potential is for baryon number 



and the baryon number of a single quark is . The diffi- 
culty in simulating this functional integral on the lattice 
stems from the fact that the functional determinant is 
not generally real and positive. From the perspective of 
QCD inequalities, however, this is not a bug, but a fea- 
ture; it allows one to bound the partition function from 
above: 



Zb{T,iib) < 
Jd[A] 
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In order for inequality |(2J) to be useful, its right-hand 
side needs to be expressed in terms of a physical quantity. 
Fortunately it can be related to the free energy density of 
QCD with an isospin chemical potential |lOj . An isospin 
chemical potential term is of the form Hiq^o^-q, which 
implies that the functional integral for the grand parti- 
tion function Zi{T,fij) = exp (—jSV^ G j(T, fij)) is given 
by 



Zi{T,ni) = 
X det(^-|-m 



d[A] e 



7o)det(^+m+^7o) (3) 



where the two functional determinants are for the two 
fiavors of quark and the opposite signs of the [ii terms 
refiect the opposite values of 1^ for the two flavors. To 
proceed we use the fact that 



75(^+m + y7o)75 = (-^+TO-y7o) 
= (^+m-^7o)t 



(4) 



where the last equality exploits the fact that in Euclidean 
space Ip is anti-Hermitian while the other two opera- 
tors are Hermitian. Exploiting the cyclic property of 
the determinant allows one to write the second func- 



tional determinant in eq. © as det(^-f 
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det (75(^-1- m -I- ^ 70)75) and using equation eq. @ 
then gives 



det(^+m+ ^70) = 

(det(^+m- ^70) 

Combining eq. ( El with eq. (j3l yields 
Zi{T,fii) = 



(5) 
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Finally, inequality together with eq. ® implies that 
Zi{T, 2^) > Zb{T, hb)- This relation together with the 
definition of the free energy requires that 

2Ais, 



GBiT,fiB) > Gi{T, 



(7) 
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which is the first principal result of this letter. 

Although this result was derived for two flavor QCD, 
the argument goes through for QCD with two degenerate 
light flavors and additional heavy flavors, provided the 
chemical potential term is understood as being the chem- 
ical potential for the up and down quark contributions to 
the baryon number rather than the full baryon chemical 
potential. The only change in the argument needed for 
this more general case is to include the functional de- 
terminant for the heavy flavors in all of the functional 
integrals. As the chemical potential does not apply to 
these heavy flavors, the functional determinants are the 
same as at /i = and hence are real and non-negative. 
Thus, they do not alter the preceding inequalities. This 
more general case is significant as in nature QCD has two 
light quarks which are nearly degenerate and additional 
heavy flavors. The result also applies in the general case 
to the full baryon chemical potential if one is in a regime 
in which the sjos = cjqc — Iryob — tjot = 0, since in this 
regime the total baryon number comes from up and down 
quarks. Such a regime occurs at zero temperature pro- 
vided the chemical potential is below the critical chemical 
potential for strangeness condensation to occur. 

A bound on the nucleon mass may be derived from 
inequality using thermodynamic arguments. The 
bound applies to QCD with two degenerate light flavors 
and any number of heavy flavors. To begin, note that in- 
equality lO holds at all temperatures, including T = 0. 
At zero temperature there are no thermal fluctuations 
and the system is in a single quantum state which min- 
imizes the free energy G — H — fj,N, where /i is the rel- 
evant chemical potential (either isospin or baryon) and 
N — VsP is the related particle number. The role of 
the chemical potential at zero temperature is simply to 
alter the relative free energies of the various quantum 
states. Thus increasing the chemical potential from zero 
at T = will have no effect until it is large enough so 
that the free energy of some other quantum state drops 
enough to equal that of the true vacuum. Accordingly, 
there is a critical value for the absolute value of the chem- 
ical potential at T = below which the density is zero 
and the free energy is that of the vacuum state (which is 
conventionally taken to be zero). The critical chemical 
potentials are thus defined as follows: 
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Inequalities iQ and ©, together imply the relation 



Mb 



> 



(9) 



Consider the quantum state of a single nucleon at rest. 
This state has energy Mjv and baryon number unity; its 
free energy is = Af^r — Mb- Clearly this is less than 
or equal to zero for /i^ > Mfq. Thus, there is at least 
one state lower than the vacuum whenever jib > Miq. 
This implies that /i^ < Mm- This last inequality would 
become an equality if the system for just above its 
critical value formed an arbitrarily low density gas of nu- 
cleons (implying a second order transition). However, 
this is not what actually happens. Based on the extrap- 
olation of finite nuclei densities and masses to infinite 
nuclear matter one has a solid empirical basis to 
conclude that the transition is first order: for /i^ just 
above pfg the system has nonzero energy and nonzero 
density. Thus /i^ = Afjv ^ B where B is the binding 
energy per nucleon of infinite nuclear matter and the in- 
equality, pLB > is not saturated. In any event, the 
inequality ^% — together with inequality © yields 
a bound on the nucleon mass. 



Mn > 



(10) 



For inequality H10|l to be useful we need to know ^j. 
For sufficiently small values of the quark mass we know 
that chiral perturbation theory accurately describes low 
energy excitations of the QCD vacuum. In chiral per- 
turbation theory, the phase transition associated with 
increasing /i/ is second order and amounts to pion con- 
densation, which implies /ij = |l6l |. More generally, 
/ij corresponds to the state in QCD with the lowest en- 
ergy per unit isospin. Let us denote such a state as X. 
There are two possibilities for X. If the transition is sec- 
ond order (as it is in chiral perturbation theory), then X 
is a single particle state; such a state is clearly at zero 
momentum so the energy can be identified as the mass of 
the particle. In this case, we can denote the mass of the 
state X, as mx, and its isospin as , Ix- The inequality 
(irU)! becomes 



Mn > 
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It is straightforward to bound the from above by 
the nucleon mass using a simple variational argument. 



Although it is strongly believed that in nature the tran- 
sition is second order, we do not require this in obtaining 
an inequality. Assume for the moment that the transi- 
tion were first order. In this case the state of minimum 
energy per isospin would in fact be infinite isospin mat- 
ter (in analogy to infinite nuclear matter), i.e. a state of 
uniform isospin density and uniform energy density fill- 
ing all space. However, if this were the case, one could 
construct single particle states by taking large chunks of 
isospin matter (in analogy to large nuclei in a world where 
the electro-magnetic interaction was shut off). The en- 
ergy of such a state will have a bulk contribution scaling 
with the volume of the chunk and equaling fJ-'jl^ where /3 
is the isospin of the state. Corrections to this bulk value 
scale as the surface area and hence go like I^^"^. Thus, by 
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making arbitrarily large but finite one can can create 
single particle states whose energy /isospin is arbitrary 
close /^J. Identifying such a state as X one finds that 
there exists an X for which can be made arbitrar- 
ily close to ^j. Accordingly, regardless of whether the 
transition is first or second order one can conclude that 
inequality H10|l implies there exists in QCD some single 
particle state X such that inequality Hll|l holds. 

Note, the state X cannot be the nucleon itself as this is 
inconsistent with the inequality above; hence, inequality 
makes a nontrivial prediction. 

A few comments on the significance of these results 
is in order. Consider the implications of inequality iQ. 
As noted previously, QCD at finite baryon chemical po- 
tential cannot be simulated on the lattice using Monte 
Carlo methods. It is generally believed that asymp- 
totically high densities can be treated analytically from 
QCD using weak coupling, but nonperturbative, tech- 
niques along the line of BCS theory 0, 0|. However, 
it is also generally believed that these densities are ex- 
traordinarily high and the regime where these methods 
work is unlikely to be relevant in laboratory experiments 
or in stellar physics. Thus, studies of phenomenological 
significance have been based on various models |^ Il2l|. 
Ideally, such models should be as constrained as possi- 
ble from QCD. Inequality provides one possible basis 
for such a constraint. While the left-hand side of the in- 
equality will be given by the model, the right-hand side 
is tractable in lattice QCD. The key point which allows 
such lattice simulations was recognized several years ago 
[lO] and is simply that the integrand in the functional 
integral for Zi in eq. |j2Jl is positive definite and hence 
is amenable to lattice studies. Lattice studies have been 
done both in the quenched approximation '14*1 and in- 
cluding dynamical quarks ilSj . It is probably sensible 
to consider the dynamical lattice studies as being rather 
preliminary as they have been done on small lattices. 
However, as the lattice calculations improve, they may 
provide a strong constraint on models of QCD at finite 
baryon chemical potential. 

Next consider inequality 1)11(1 . It is very similar to the 
bounds derived by Nussinov for a wide class of QCD- 
inspired potential models and to the bound derived 
by Nussinov and Sathiapalan in large Nc QCD 6]. The 
present result suffers in comparison to these results in 
that instead of bounding the nucleon mass by Ncm-,^12 
the bound is with respect to M^v > where X is 

the state in QCD with lowest mass per unit isospin. In 
practice this disadvantage is fairly small — it so happens 
that in QCD the state with the lowest mass per unit 
isospin is the pion. The present bound has a clear ad- 
vantage over the the results in refs. 2, 6] in that it holds 
for QCD itself and is not restricted to the large Nc limit. 
Similarly, the present bound is much stronger than Wein- 
garten's original bound 3] in that it holds for QCD with 



two degenerate light flavors and does not require the ad 
hoc assumptions about the quark propagator in the pres- 
ence of a background field. In any event the inequality 
((ll|l is true phenomenologically. There does exist a state 
in QCD for which Mn > (3Mjf)/(2/x) namely the case 
where X is the pion. Overall, one should view any bound 
on the nucleon mass in terms of other physical observ- 
ables directly from a first principles treatment of QCD 
as being highly non trivial. 
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